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Abstract 

The interaction of two identical three-level atoms of the types V, H and A with a quantized 
cavity field as well as a driving external classical field is studied. Under two certain unitary 
transformations, the system is converted to a typical form of the Jaynes-Cummings model for 
two three-level atoms. The exact analytical solutions of the wave function for different considered 
atom-field systems are exactly obtained with the help of the Laplace transform technique, when 
the atoms are initially prepared in the topmost excited state and the quantized field is in a coherent 
state. In order to examine the nonclassicality features of the deduced states, the dynamics of the 
entanglement between subsystems is discussed via two well-known measures, namely, von Neumann 
entropy of the reduced state and negativity. In addition, we pay attention to the temporal behaviour 
of quantum statistics of the photons of the field and squeezing phenomenon. Meanwhile, the 
influence of the external classical field on the latter physical quantities is analyzed in detail. The 
results show that the mentioned quantities can be sensitively controlled via the external classical 
field. Also, numerical computations imply the fact that the nonclassicality features in H-type three- 
level atomic system is more visible than the other two configurations. In addition, it is shown that 
in the particular case of A-type atomic system, the rank of the reduced density matrix of the 
three-level atoms is no larger than three, so that negativity fully captures the entanglement of this 
system and that such entanglement is distillable. 
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I. INTRODUCTION 


The light-matter interaction is an essential concern in optical physics. A simple paradigm 
of this interaction contains a two-level atom conpled to a single-mode qnantized radiation 
held in an optical cavity. Whenever the strength of the atom-held coupling is far smaller 
than the held frequency, the rotating wave approximation (RWA) is aj^licable and the sys¬ 
tem is described by the well-known Jaynes-Cummings model (JCM) [l|, 1^. The simplicity 
of this model and its potential applications for more complicated and generalized atom-held 
systems, together with revealing some non-trivial phenomena and extraordinary character¬ 
istics, such as the collapses and revivals of Rabi oscillations in atomic population inversion 
[^, placed the JCM at the heart of this area of research in quantum optics. Experimentally, 
this model can be realized when atoms are coupled to a nanomechanical oscillator jj] and 
nuclear spins interacting with a magnetic held Many generalizations of the JCM have 
been proposed in various ways, for instance, considering diherent initial conditions jb], enter¬ 
ing the ehects of dissipation and damping in the model j^, considering intensity-dependen 


coupling 


11| . adop 


and a multi-atom 


16 


ing a multi-level atom 12l-ll4| as well as a multi-photon transition 151 


17 |. 


From another perspective of this area of research, in recent decades, a lot of attention has 
been paid to the study of diherent conhgurations of three-level atoms (U, S and A type) 
. For example, there exist many theoretical works containing the interaction between a 
three-level atom and a single-mode cavity held 18l-l21|. In some of them, some nonclassi- 


cality features have been examined by considering the resonance condition between atomic 
transitions and quantized held frequency. The nonclassical properties of three-level atomic 
systems have been well investigated in order to understand the^antum coherence phenom¬ 
ena such as electromagnetically induced transparency (FIT) 2^, lasing without inversion 


23l |. and coherent trapping [2^. In detail, the U-type three-level atoms are extensively 


used in studying nonclassica’ 


ity features such as quantum beats 


24l |. quantum Zeno ehect 


25| and quantum jumps |26|. The application of S-type atoms is well established in the 


coherent population trappin 
laser cooling in trapped ions 


27l | and also the experiments which are designed to achieve 
Also, the A-type atomic systems have been widely utilized 



in representing the coherent phenomena such as FIT 22j and stimulated Raman adiabatic 


passage 


29]. 
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In particular and in direct relation to the present work, the dynamics of a physical system 
governed by the JCM can lead to the highly atom-field entangled state, almost as though the 
atom and the held can form some sort of “molecule” [^; the case that is obviously revealed 


through the vacuum Rabi splitting 


3l| . It may be noted that, although, such experiment 


regarding the observation of a molecule is not usually a routine work, however, this may be 
naturally probed through an external held. Furthermore, strong driving helds have shown 
other attractive ap plic ations in the atom-held interaction such as detection of Fock states of 
the rac.: 


states 


la 


34 


ion held 


32| . quantum-phase gate j33|, generation of the multi-partite entangled 


36l | and controlling nonclassical properties of system 37|. Henceforth, the driving 


JCM (DJCM), that is, the JCM which contains an external driving held, has attracted much 


attention in recent decades 


ss-dil . 


The model for a two-level atom interacting with exter¬ 


nal quantum and classical electric helds through a parametric frequency converter has been 
presented in |^. In this attempt, the system is detracted to an ehective JCM by adequate 
adjustment of the held coupling in the frequency converter. In addition, it has been shown 
that the dynamics of physical quantities related to the atom and the atom-photon entan¬ 
glement can be controlled by the classical held. The implementation of a strongly driven 
one-atom laser, based on the oh-resonant interaction of a A-type three-level atom with a 
single-mode cavity held and three laser helds have been schematically proposed in The 
authors showed that the system can be equivalently described by a two-level atom reso¬ 
nantly coupled to the cavity held assisted by a strong ehective coherent held. Altogether, 
in all above discussions, it seems that the driving external held is a suitable parameter for 
controlling the nonclassicality features, especially entanglement of the atom-held systems. 
In this paper, we aim to study the problem of two identical three-level atoms in three dif¬ 
ferent types (R, S and A conhgurations) interacting with a single-mode cavity held in the 
presence of an external classical held through which new classes of entangled state can be 
generated. In working with the mentioned conhgurations, we are able to transform the 
interaction formalism to the generalized JCM by applying two appropriate unitary trans¬ 
formations. The explicit form of the entangled state vector of the whole system can then 
be exactly obtained by using the time-dependent Schrodinger equation with the help of the 
Laplace transform technique. In this respect, we take the atoms to be prepared in their 
higher excited states and the held is supposed to be in a coherent state. Then, briehy the 
main goal of this paper is to investigate the ehect of external classical held on the entan- 
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glement dynamics between subsystems and some of the well-known nonclassicality features. 
To achieve this purpose, the degree of entanglement (DEM) between subsystems through 
von Neumann entropy of the reduced state (to investigate the atom-held entanglement) and 
negativity (to obtain the atom-atom entanglement) are numerically studied. Also, the time 
evolution of quantum statistics and squeezing are examined in detail. We show that the 
dynamics of the mentioned physical quantities can be tuned by the strength of the external 
classical held. 

In order to give more explanation about our motivations, it is instructive to give a few 
words regarding the signihcance and the notability of the considered systems containing the 


three-level atoms. In quantum information processing (QIPVt 


outstanding advantages in comparison with two-level ones 


le three-level systems possess 


44| . In this regard, the optimal 


eavesdropping in quantum cryptography with three-dimensional systems has been studied in 


45| in which the authors found that the three-dimensional scheme ohers higher security than 


the two-dimensional systems. Kaszlikowski et al investigated the general case of two entan¬ 
gled quantum systems dehned in d-dimensional Hilbert spaces, or ‘qudits’ and have shown 
that violations of local realism are stronger for two maximally entangled qudits (3 < d < 9) 


than for two qubits 


46| . Also, based on the Greenberger-Horne-Zeilinger (GHZ) theorem. 


the conflict of local realism and quantum mechanics for three or more qubits has been re¬ 


ported much sharper than for two qubits 


47( 1 . In this case, there exist some suggestions 


which indicate the reduction of this conflict due to increasing d (dimension of Hilbert space) 


48 


50|. So, it seems that d-level quantum systems, or qudits, may be better candidates 


to be utilized in the theoretical/experimental observations. The mainspring of this usage 
is to increase the available Hilbert space with the same amount of physical resources 51 |. 
Accordingly, three-level atoms have received noticeable attention in the studies that concern 
with the atom-held interaction 


52 


55| . In particular and in direct relation to the consid¬ 


ered model in this paper, quantum information processing using superconducting qubits 
has made outstanding advances in the past few years [56|, l57(]. One-qubit and two-qubit 
quantum circuits have been realized experimentally in superconducting systems. One of the 
most important issues in quantum information processing is how to couple two qubits, which 
has been widely studied theoretically and experimentally in superconducting quantum cir¬ 
cuits. Theoretical proposals have been put forward to selecti vely couple any pair of qubits 
through a common data bus (LG circuit or a cavity held) 58|, l59j. Liu et al. experimentally 
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presented a proposal to achieve a controllable interaction between flux qubit by virtue of 


time-dependent electromagnetic field 


60| . Also, a coupling (decoupling) method between a 


superconducting qubit and a data bus that uses a controllable time-dependent electromag¬ 
netic held (TDEF) have been studied theoretically in This reference indicates that, by 
choosing appropriate parameters for the TDEF, the dressed qubit (qubit plus the electro¬ 
magnetic held) can be coupled to the data bus and, thus, the qubit and the data bus can 
exchange information with the assistance of the TDEF. The superconducting qubit circuits 
generalized to superconducting qudit circuits for more complex quantum computational ar¬ 


chitectures, and for richer simulations of quantum mechanical systems 


62| . The considered 


model in this paper may be supposed as a superconductor with two qutrits (three-level 
atoms), LC circuit (single-mode cavity held) and TDEF (the classical held). In this case, 
the entanglement between two qutrits as well as between two qutrits and a single-mode held 
can be controlled by external classical held. 

The reminder of the paper is structured as follows. In the next section, the model containing 
all existing interactions is introduced and then by applying two distinct unitary transfor¬ 
mations, the model is reduced to the typical form of the generalized JCM. In section HTTl 
the state vector of the whole atom-field systems is analytically obtained. Section IIVI deals 
with examining the effect of external classical held on the DEM between subsystems via von 
Neumann reduced entropy and negativity. In addition, in order to study the nonclassicality 
features of the obtained states, Mandel parameter and quadrature squeezing are respectively 
investigated in sections IVl and IVTl Finally, the main results of the paper are summarized in 
section IVIIl 


II. DESCRIPTION OF THE MODEL 


We consider two identical three-level atoms (labeled by A and B) in three different 
configurations (namely V, S and A types as depicted in figure 1), with states |1), |2), and |3) 
and their corresponding energies Ui, U 2 , and The atomic system is driven by an external 
classical held with frequency cUc, and is coupled to a single-mode quantized radiation field 
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with frequency Uf. Accordingly, the total Hamiltonian can be appropriately described by 


3 

H= ^ + cu/a'^'a + ^ 

j=A^B i=l j=A,B 

+ 92 

j=A^B j=A,B 

+ X2YI ( 1 ) 

j=A,B 

where an = |z)(i| is the atomic projection operator, a and d'*' are respectively the bosonic 
annihilation and creation operators of the held, gi and Aj, i = 1,2 represent the coupling 
constants of the interaction of the atoms with the quantized radiation and with the classical 
driving helds, respectively. Also, the values of (di, d 2 ) for the three conhgurations are given 
by (di 3 ,d 23 )y, (di 2 ,d 23 )H and (di 2 ,di 3 )A. In the rotating reference frame with frequency 
ojc and under a unitary transformation Ui{t) = exp[—iu)ct{aXa + J2j=AB (*;<'■'+*;«>))]. the 
above Hamiltonian can be transformed to 


H, = Ul{t)HU,{t) - 

at 

= Ai d/ + A2 a^ + Ad^d 

j=A,B j=A^B 

+ 9i 51 (hP® + +92 51 (+ 2 ’'“+ + 2 ™“’) 

+A/x;(hf'+9r)+A2 5i(9“+d“). 

j=A,B j=A,B 


( 2 ) 


where the values (d^,d 2 ) of the three-level atom are in the form (dii,d22), (dn, —d 33 ) and 
(—d 22 ,—d 33 ), respectively for V-, H- and A-type atoms. The parameter A = Uf — Uc is 
the detuning parameter between the cavity held and the classical driving held. Also, the 
detuning parameters between the classical held and the atoms are given by Ai = (cji — Hi) — 
oJc and A 2 = (H 2 — a; 3 ) — Uc, in which the values of (Hi, H 2 ) in three conhgurations read as 
(a; 3 ,a; 2 )y, (<+.; 2 ,<^ 2 )h and (a; 2 ,n;i)A. It is seen that by applying the mentioned transformation, 
the time-dependent exponential terms in the Hamiltonian ([T]) have been clearly eliminated. 
Now, for simplicity, let us follow the problem in the resonance conditions. Also, without 
loss of generality, we suppose Qi = 92 = 9 and Ai = A 2 = A. Under these assumptions, we 
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can recast the Hamiltonian (|2]) as follows 


Hi = g + a2^a + a) + a)) 

j=A,B 

j=A,B 


( 3 ) 


If A = 0, then Eq. ([3]) describes the generalized JCM for the interaction between two three- 
level atoms and a single-mode quantized cavity field, while the case A 7 ^ 0 corresponds to 
generalized DJCM for two three-level atoms. Now, in order to analyze the dynamics of the 
considered system with the Hamiltonian ([3]), we utilize the probability amplitude method. 
Altogether, it is not still an easy work to solve the above system analytically. This is due 
to the existence of the external classical field in addition to the terms which are connected 
to the interaction between two three-level atoms and single-mode quantized field. Therefor, 
before using this approach, we try to reduce the Hamiltonian ([3]) to the typical form of the 
generalized JCM for two three-level atoms. This goal will be achieved by introducing the 
following displacement operator 


D{'y) = exp(7a^ — 7*0), 7 

which satishes the identity 

D{'y)aD\'y) = 0 — 7. 


A 

9 


( 4 ) 

( 5 ) 


In this case, by applying the unitary operator given in (j4]), the Hamiltonian ([3]) is converted 
to the Hamiltonian of the form 


H 2 = D{'y)HiD\'y) = g ^ df'^^d'*' -h (6) 

j=A,B 

Briefly, up to now we used the two transformations (the local transformation tJi{t) and the 
displacement operator D{'y)) for reducing the Hamiltonian ([T]) to the typical generalization 
of the JCM for two three-level atoms. It is worth noticing that, the physical properties as 
well as the initial conditions of any system are preserved under the local transformation. 
But the second transformation (displacement operator) changes the initial condition and 
the physical features of the system. Henceforth, by using the Hamiltonian ([H]) and defining 
the state vector |' 02 (^)) = one may consider the time-dependent Schrodinger 

equation corresponding to \'ijj 2 {t)) as id\'ijj2{t)) / dt = H 2 \'^ 2 {t)), from which one easily obtain 
I'lpit)) = {■j)\'ip 2 (t)). In the next section, we are going to evaluate the state vector \'ip 2 (t)) 

for different configurations of the three-level atoms. 


III. THE ENTIRE STATE VECTOR OF THE SYSTEM 


The main goal of this section is to obtain the state vector of the considered systems. 
Before achieving this purpose, it is necessary to determine the initial conditions of the 
atoms as well as the held. We consider the held to be initially in the coherent state and 
suppose that the atoms enter to the cavity in the upper exited state, i.e.. 


|V’(t = 0)) 




a) = exp 



(7) 


Hence, the state vector = 0)) may be obtained by the following relation 


\Mt = 0)) = I)(7)|l, l)l«) = |1,1)^(7)^(«)|0) = |1,1)1/3), (8) 


where we have used the identity 

D{'y)D{a) = exp{i^{'ya*)) D{I3), (9) 


so that /3 = a + 7 . The factor exp{i‘^{'ya*)) is dropped because the parameters 7 and a 
have been supposed to be real. By considering this initial condition for \'ijj 2 {t)), the wave 
function for the considered systems are obtained separately, in the next subsections. 


A. Two V-type three-level atoms 


The wave function \'ijj 2 {t)) at any time t for two H-type three-level atoms may be written 

00 / 

\i^ 2 {t))v = (Ci{n,t)\l,l,n) + C2{n,t){\l,2,n) + |2,l,n)) 

n=0 ^ 

+ + l,t)(^|l,3,?7- + 1) + |3,1,?7. + 1)) 

+ C*4(?7- + l,t)(|2,3,?7. + 1) + |3,2,?7. + 1)) 

+ C 5 (n, t) |2, 2, n) + Ce{n + 2, t) |3, 3, n + 2)^ , (10) 


where the coefficients Ci,i = 1 , 2 . .. 6 , are the unknown probability amplitudes that should 
be determined. By inserting the assumed wave function OB into the time-dependent 
Schrodinger equation together with considering the Hamiltonian ([ 6 ]), the following coupled 
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differential equations for the probability amplitudes may be found: 


dCi{n, t) 
dt 

dC2{n,t) 

dt 

dC^in + l,t) 

dt 

dCi{n + 1, t) 
dt 

dC^{n, t) 
dt 

dCg[n + 2, t) 
dt 


= -iViC‘i{n + l,t), 


— —iVi I C^iji + 1, t) + + 1, t) I, 


= -iVi Ci{n,t) + C 2 {n,t) - iV 2 C^{n + 2,t), 


= -iVi C2{n, t) + C'5(n, t) - iV2C^{n + 2, t), 


= —2iViCi{n + 1, t), 


— —2*1/2 I C^iji + 1, t) + + 1, t) I, 


( 11 ) 


where Vi = V{n + 1) and 1^ = l^(n. + 2) and V (n) = g^/n. After some lengthy calculations, 
we obtain the probability amplitudes via the Laplace transform techniques as below: 


C,(., 0 ^ 0 ) ( 2 cos(^^,) + 1 . 




(12a) 


C 2 {n,t) = --C'i(n,0)- — - 


(12b) 


C 3 (n + l,t) = ^C'i(n, 0)| \/2sin(\/21/it) + 


(12c) 


C 4 (?t, + 1, t) = -Ci{n, 0) I \/2sin(\/2Vit) 


Vi sm{2'dt) 


(12d) 


C,{n,t) = lc,(n.O) ( - 2cos(V21/.() + (i2e) 


Cq^u + 2,t) = -Ci{n, 0) 


1/11/2 sin^('dt) 
^2 ^ 


(12f) 


where id = \/Vx + ^^*1 C'i(^, 0) = exp(—|/9p/2)/3”'/\/rd determines the probability of 

the initial field state. 
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B. Two H-type three-level atoms 


In a similar manner, the state vector |t/’ 2 (^)) for two S-type three-level atoms is given by: 

OO / 

|' 02 (t))H = ( Ci{n,t)\l,l,n) + C 2 {n + l,t){\l,2,n + 1) + |2,l,n-M)) 

n=0 \ 

-|- C'3(?7,-|-2,t)(|l,3,?T.-|-2) -1- |3,1,?t.-|-2)) 

-|- 04(71 -|- 3 , t )(| 2 , 3 ,? 7 .- l - 3 ) | 3 , 2 ,? 7 .-|- 3 )) 

-|-C*5 (?t.-|-2, t) |2, 2,77.-|-2)-|-Cg(?7.-|-4, t)|3, 3, n-|-4) 1. (13) 


Following the same procedure as subsection IIII A 1 the values of the coefficients C* are ob¬ 
tained as below: 

C'i(n, t) = ( (x2 - X4)r] (2F^^X2 - (51x4) cos(/ 3 it) 

X2ri 


( 2 V 4 X 2 - (51x4) cos(/52t) , 


( 14 a) 


02(71 -|- 1 , t) — 


iOi(n, 0 ) 


(x4 - 2^^V^ )^2sm(/3it) 


(x4 - 2/32foi )/9i sin(/32t) , 


( 14 b) 


03(71 + 2,t) = I -X57]- ((5^X5 - V1V2X2) cos(/ 3 it) 

X27] 


+ (/3lx5 - V1V2X2) cos(/32t) , 


^4(77 + 3 ,t) = I sin(^it)/^i - sm(( 52 t )/(32 ), 


( 14 c) 

( 14 d) 


6*5(77 + 2 ,t) = I -X37]- ((31x3 - V1V2X2) cos(( 3 it) 

X27] 


+ ((31x3 - V1V2X2) cos((32t) , 


6*0(77 -|- 4 , t) — 


0:16*1(77,0), 2 


X27] 


7]- ( 3 f cos(/52t) + /da cos(( 3 it) , 


(foe) 

(fof) 
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with 


xi = 6i/ii/2^3^4, a;2 = + 6V2%^ 

X 3 = 2(i/2 + + 3 (Vi + If), 0:4 = + 4 V,%^, 

X5 = 2V1V2V4, r] = Jxl - 4x2, 


ft = 


ft = 1 /^. 


Vj^V(n+j), j = l,2,3,4, V{n)=g</n,. 


(16) 


C. Two A-type three-level atoms 


Finally, the wave function of a system containing two A-type three-level atoms and a 
single-mode cavity field can be written as follows 


\Mt))A = Y1 


Ci(n, f)|l,l,n) + C 2 (n + l,f)(|l,2,n + 1) + |2,1,77,+ 1)) 

+ C^iji + l,f)(|l,3,?7. + 1) + |3,1,?7- + 1)) 

+ C' 4 (n + 2,f)(|2,3,n + 2) + |3,2,n + 2)) 

+ C*5(?T. + 2,t)|2,2,?7, + 2) + Cq^ti + 2,f)|3,3,?T. + 2)) 


(16) 


where Cj are the time-dependent probability amplitudes which must be obtained. Similarly, 
we arrive at 


Ci{n,t) = Ci(n,0) 


Vi + Vi cos( 21 / 3 t) 


1/2 

*^3 


(17a) 


C 2 {n + l,t) = C 3 {n + l,t) = -iCi{n,0) 


Vi sin( 21 / 3 f) 
2 K ^ 


(17b) 


C*4(?t. + 2, f) — C^iji + 2, f) — Cq{ji + 2, t) 

1/11/2 sinfl/gt) 


= -Ci(n,0)- 

where Vi and V 2 have been previously defined in ([15]). 
Here, it ought to be mentioned that, the relation 

ji^(t)}v,E,A = ^(- 7 )lV' 2 (t))ys,A, 


^3^ 


_ A 
^ 9 


(17c) 


(18) 


12 












allows us to write the solution of DJCM (Eq. ([T])) for two V-, S- and A-type three-level 
atoms, explicitly. For example for two A-type atoms the state vector is in the following form: 


IV'W)a = 


Ci{n,t)\l,l)\ - 7 ;n) 


-|- C 2 {n + 1, f) ^|1, 2)1 — 7 ; n -|- 1) -|- 12 ,1) I — 7 ; n -|- 1) j 
+ C 3 (n + - 7 ;n + 1) + |3,1)| - 7 ;n-f- l)j 

-|- C 4(^71 -|-2,f)^|2,3)| — 7)^“I“2) -|- |3,2)| — 7j'^“l“2)^ 
+ C' 5 (n + 2 ,f)| 2 , 2 )|- 7 ;n + 2 ) 


+ C' 6 (^ + 2, f) |3, 3) I — 7 ; n -f 2)) 


(19) 


where \ —^]n + j), j = 0,1,2, are the displaced number states. Similar expressions can 
be simply obtained for other two types of atoms. Therefore, our proposed model can also 


be considered as a novel scheme for the generation of displaced number states 


By 


considering similar approach for two other cases {V- and S-type three-level atoms), the 
displaced number states can be generated, too. Anyway, we are now ready to study the 
nonclassical properties of three different types of the atom-field system by emphasizing on 
their entanglement properties. For this purpose we will pay attention to entanglement, 
photon statistics and quadrature field squeezing. 


IV. ENTANGLEMENT 

Entanglement is the noticeable feature of quantum states which demonstrates correlations 
that cannot be classically accounted. The first investigations of entanglement date back to 
1935, focusing on surprising consequences of the quantum description of nature [GJ, l65|. 


Entangled qubits are an urgent resource in many quantum information applications such as 
quantum computation and quantum communication [ 66 | , quantum metrology |67| , quantum 
cryptog raphy | 68 l |. quantum teleportation |69| and other applications in quantum technology 
70|, l7l| . Recently, much attention has been paid to the generation of quantum entangled 
states. A well-known source for the generation of such states is the atom-held interaction 
process, using different models of interaction. Accordingly, the most interesting aspects of 
the JCM and its generalizations, which has received much attention, is the possible existing 


of entanglement between different subsystems 


72 


74| . So, we now pay attention to the 
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evaluation of the entanglement dynamics of the obtained states. To achieve this goal, several 
suitable measures of DEM such as von Neumann entropy [^, entanglement of formation 


0 . 


concurrence 


76 


77| and negativity 


78| have been proposed. In this section we apply 


von Neumann entropy and negativity to investigate entanglement dynamics of atom-held 
and atom-atom, respectively. It ought to be mentioned that, while the von Neumann entropy 
is a good measure for the atom-held entanglement, this measure is not appropriate for the 
calculation of the DEM between the two atoms (this is due to the fact that in this case, 
the system (the two atoms) is a mixed state). Also, the evaluation of DEM between “two 
atoms” and “held” by negativity is not an easy task (its complexity arises from the fact that 
to achieve this purpose we are left with the inhnite dimensional Hilbert space related to the 
helds). 


A. von Neumann entropy 


We have assumed that two three-level atoms and the coherent held are initially in a pure 
state. So, the considered systems can be regarded as the bipartite systems, consisting of two 
three-level atoms as the hrst subsystem and the radiation held as the second subsystem. For 
such systems, the von Neumann entropy is a suitable measure to obtain the DEM between 
subsystems 79|. Quantum mechanically, the von Neumann entropy for a quantum system 
with the density operator p is dehned as S' = —Tr(plnp). If p describes a pure state, 
then S' = 0, and if it represents a mixed state then S' 7^ 0. For the considered atom-held 
systems, the entropy of the held or equivalently the atoms is a good measure to realize the 
amount of entanglement; higher (lower) entropy means the greater (smaller) DEM. Before 
obtaining the reduced entropy of the held and the atom, it is worth to pay attention to 
the important theorem of Araki and Leib 8^ . According to this theorem, for any bipartite 
quantum system, the system and subsystem entropies at any time t are bounded by the 
triangle inequality |S'yi(t) — S'i?(t)| < S < |S'A(t) + S'i7’(t)|, where S'a and Sp represent the 
entropies of the atom and held, respectively, and the total entropy of the atom-held system 
is denoted by S'. One immediate consequence of this inequality is that, if at the initial time 
the whole system is prepared in a pure state, the total entropy of the system is zero and 
remains constant when the involved subsystems are isolated from their environment. This 
implies that, if the system is initially in the pure state. Sap = 0, at any time t > 0 the 
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field and atomic entropies are equal Ml- Therefore, instead of the calculation of the field 
entropy, we can evaluate the entropy of the atoms. According to the von Neumann entropy, 
the entropies of the atom and field, when treated as a separate system, are dehned through 
the corresponding reduced density operators as 


SA{F){t) = —TTA(F){pA{F){t) \n pA{F){t)), 


( 20 ) 


with pA{F)(t) = Tri?(A)(|'^(t))(t/’(t)|), as the reduced density operator of the atoms (held). 
Now, we turn our attention to discuss the DEM between two atoms and quantized held for 
three diherent systems through the von Neumann entropy. We assume that the considering 
systems start from a pure state, so the entropy of the held/atom may be expressed by the 
following relation 

9 

DEM(t) = Spit) = SA{t) = (21) 

i=l 

where denote the eigenvalues of the reduced density matrix of the atoms, which can be 
obtained numerically for V- and H-type, but analytically for A-type conhguration. Indeed, 
the density matrix of the two-atom system in the A-type conhguration has rank no larger 
than three with the associated eigenvalues evaluated, analytically, by the Cardano’s method 
as 3 

= ~Qi + \\l Qi- 3p2 cos - l)7r^ , i = 1,2,3, 


A = 0, j= 4,5, 6, 7,8,9, 


( 22 ) 


with 


and 


VJ = - cos 
3 


-1 


9piP2 — 2pf — 27p3 
2(£>;-3&)V2 


(23) 


^?l — —pll — 4(p22 + ^ 44)5 


(24a) 


Q2 — —4:{pl2P21 + Pupil + 4 ^ 24 ^ 42 ) + 4pii(p22 + P 44 ) + 16^22^44, 


(24b) 


P3 — 16pi4(p22P41 ~ P21P42) + 16pi2(P2lP44 ~ p2ipil) 
+ 16pil(p24P42 ~ P22pA4)- 


(24c) 
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Moreover, the matrix elements of the atomic density operator are as follows: 

+ CX) 

= hi = 1,2, •••6. (25) 

n=0 

where Ci{n,t) have been derived in [T71 Figure 2 shows the evolution of the held entropy 
against the scaled time gt for initial mean photon number hxed at |ap = 25 and two atoms 
prepared initially in the higher exited state. Frames 1, 2, and 3 respectively concern with the 
two V-, S- and A-types three-level atoms. Also, panels (a) in this hgure are plotted in the 
absence of classical held (7 = 0) and panels (b) and (c) show the ehect of the driving external 
held on the behaviour of von Neumann entropy (7 = 2 and 7 = 6). It can be obviously seen 
from this hgure that for three diherent conhguration of three-level atoms, the held entropy 
gets the maximum value of entanglement after the onset of the interaction. In addition, 
it can be observed from frame 1 that, entering the classical held together with increasing 
the related parameter may lead to the shift of the maxima amounts of the DEM when the 
time proceeds. Comparing the panels 1 and 2 indicates that the temporal behaviour of two 
S-type three-level atoms is qualitatively the same as two D-type ones. For the case that 
we deal with A-type three-level atoms, by considering the ehect of external classical held, 
it seems that the temporal behaviour of the DEM behaves oscillatory specially in the case 
7 = 6. Finally, by focusing on the ehect of external classical held on the behaviour of the 
von Neumann entropy, it is found that the existence of the classical held may increase the 
maximum values of DEM between two atoms and held in three considered systems. 


B. Negativity 


In the present subsection, we apply the negativity measure for the investigation of the 
DEM between two atoms. Among all entanglement measures, negativity surely is the best 
known and most popular instrumentation to specify bipartite quantum correlations 83|. 


It is easily evaluated for arbitrary states of a composite system and so can be applied 
to study entanglement in many diherent situations. Historically, this qua ntity , which is 


based on the Peres-Horodec 
used by Zyczkowski et al 


d criterion for the separability of a state 


m 


85|, was hrst 


and subsequently introduced by Vidal and Werner as a new 


entanglement measure 78(]. The positive partial transpose (PPT) is a necessary and sufficient 


condition for separability of 2 x 2 and 2 x 3-dimensional mixed states 


M 


, arbitrary dxd' 
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pure state, all Gaussian states of 1 x n mode continuous variable systems [87|, 


and also 


any state p which is supported on d x d' systems (d < d') and with rank r(p) < d' |89|, but 
it is only sufficient, in general, for other systems. The negativity for a bipartite quantum 
system with d x d' {d < d')-dimensional Hilbert space Ha described by the density 

matrix p, is given by (d < d') 


U{p) = 


^Tb I 


d — 1 


(26) 


where p^® is the partial transpose of the state p with respect to subsystem B and ||M||i = 
TiVmIM is the trace class norm of the operator M, which reduces to the sum of the 
absolute value of the eigenvalues of M, when M is Hermitian. The matrix p is a positive 
operator with trace one, i.e. Tr(p) = 1. Also, for the partial transpose of this matrix we have 
Tr(p^s) — Since the partial transpose of density operator might have the negative 

eigenvalues, the trace norm of p^® can be written in the following form 


11/5^* 11. = E l'“l = Z - 2 E = 1 - 2 E (27) 

i i i i 

where pi and correspond to the positive and negative eigenvalues of p^®, respectively. 
Here, we are going to evaluate the DEM between two three-level atoms ((3 x 3)-dimensional 
Hilbert space) in the considered systems. 

In Fig. 3, we have plotted the time evolution of the negativity as a function of the dimen¬ 
sionless time gt for the same parameters as in Fig. 2. One can see from this hgure that the 
temporal behavior of the negativity in various conditions represents irregular oscillations 
between minima and maxima values. Frame 1 of this figure shows the DEM between two 
H-type three-level atoms with and without the external classical held. According to these 
hgures, it may be seen that by increasing the amplitude of the external classical held, the 
maxima values of negativity are revealed with time delay. For two S-type three-level atoms, 
the outstanding ehect of the classical held is decrement the sustainment time of the maxima 
values of the negativity. The same behavior can be seen in frame 3 for two A-type three- 
level atoms. Finally, by considering the presented results depicted in Fig. 3, it is found that 
the DEM between S-type three-level atoms in the presence or absence of the classical held 
is greater than two other conhgurations. Also, for diherent types of three-level atoms the 
negativity can be managed by the driving external classical held. On the other hand, the 
addition of the classical held causes a slight increase in the maximum values of the DEM 
between two atoms. 
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V. PHOTON STATISTICS: THE MANDEL PARAMETER 


Sub-Poissonian statistics is a striking feature of nonclassical states. To investigate the 


statistical properties of any system the Mandel parameter is a suitable measure 
parameter has been dehned as the following form 


90[. This 


Q{t) = 


(n^) — {nY 
(n) 


1 . 


(28) 


This quantity is positive, zero and negative when the statistics is super-Poissonian, Poisso- 
nian and sub-Poissonian, respectively. The sub-Poissonian statistics is a sufficient but not 
necessary condition for nonclassicality of the held 9l|. For our considered systems we have: 


(n) = (V'(t)la^alV'(f)) = (V'2(t)|^(7)a^aT>^(7)|V’2(t)) 

= - 7)(a - 7)IV’2(i)) = 7 ^ + (V'2(t)|a^a|V’2(t)) 

- (29) 


(n^) = = (n) + {'ip2{t)\{a^ - if {a - if\'ip2{t)). (30) 

To visualize the effect of the external classical held on the photon statistics for the diherent 
three conhgurations, we have plotted the Mandel parameter against the scaled time gt for 
diherent values of the coupling parameter ratio 7 in Fig. 4. In the absence of classical held, 
the oscillations of the Mandel parameter in three conhgurations show collapse and revival 
phenomena as would be expected. Also, in these cases the Mandel parameter varies between 
positive and negative values, which means that the photons display super- or sub-Poissonian 
statistics for diherent intervals of times, alternatively. Due to the presence of external held, 
we observe that the Mandel parameter possesses a periodic behavior in the positive region 
at the most of times. This means that in these times, the entire atom-held state of the 
considered systems has a super-Poissonian statistics. Altogether, in the presence of classical 
held and only in the beginnings time of the interaction the held has sub-Poissonian statistics. 
Increasing the value of the classical held coupling parameter leads to the decrease in the time 
average of the sub-Poissonian statistics. Finally, it is to be noted that the negative value of 
Mandel parameter in S-type is greater than its counterparts in both V- and A-types. 
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VI. NORMAL SQUEEZING OF THE FIELD 


In the present section we focus on the squeezing phenomenon. This parameter is described 
by decreasing the quantum fluctuations in one of the field quadratures with the price of an 
increase in the corresponding conjugate quadrature, such that the Heisenberg uncertainty 
principle is not violated. The squeezed light has attractive applications in optical commu¬ 


nication networks and gravitational wave detection 


92 


93| . The squeezing parameter has 


been dehned in various ways such as normal squee zing , amplitude-squared squeezing j94| . 


higher-order squeezing 


95| and principal squeezing 9g . However, in the present section we 


pay attention to the normal squeezing. 

To investigate the normal squeezing of the field, we introduce two quadrature field oper¬ 
ators X = {a + a^)/2 and y = [a — a))/2i. The system then would possess squeezing if 
one of the quadratures convinces the inequality ((Af)^) < 0.25 or ((Ay)^) < 0.25 where 
{(Afj)^) = (fj^) — Xi = X and y. Equivalently, if we define Sx = 4((Aa;)^) — 1 and 
Sy = 4:{{AyY) — 1, squeezing occurs in x (y) component if —1 < < 0 (—1 < Sy < 0). 

These parameters can be rewritten as follows: 

S', = 2{d^a) + 2^{d^) - 4(3?(a))2, (31) 

Sy = 2{d^d) - 25R(a^) - 4($>(a))^ (32) 

where (a^a) has been previously dehned in (l2^ and {a^) can be obtained in the following 
form 


(a”*) = . (33) 

Our results presented in Fig. 5 indicate the time evolution of the quadrature squeezing 
parameters S', and Sy against the scaled time gt for different types of three-level atoms as 
well as various values of the coupling parameter ratio. We can see from this hgure that 
for all cases, the squeezing exists in the x quadrature and no squeezing is occurred in the 
y quadrature. Meanwhile, squeezing arises in x quadrature only at the beginning of the 
occurrence of the atom-held interaction. Furthermore, a comparison between frame (a) and 
frames (b,c) of diherent conhgurations of three-level atoms shows that the amounts of this 
nonclassical ehect in x diminishes by increasing 7 . Also, our conclusions represent that, 
the squeezing in x component in the S-conhguration is stronger than for the other two 
conhgurations. 
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VII. CONCLUSION 


Due to the importance of the three-level atoms as well as the presence of the driving 
classical held in the atom-held interactions, in this paper, we have outlined two identical 
three-level atoms (in V, S and A conhgurations) interacting with a quantized single-mode 
held assisted by an external classical held. It is shown that the introduced system can 
be transformed to the usual form of the generalized JCM by using two appropriate unitary 
transformations. Therefore, we could solve the dynamical problem and hnd the explicit form 
of the entangled state vector of the three diherent considered atom-held systems analytically, 
by considering the atoms initially in the higher exited state and the quantized held in the 
coherent state. Next, at hrst, the quantum entanglement between the subsystems of the 
generated states are computationally evaluated by using the approach of von Neumann 
entropy (to study the DEM between two atoms and quantized held) and negativity (to 
investigate the DEM between two atoms). Then, the quantum statistics and quadrature 
squeezing of the obtained states have been numerically investigated. In each case, we studied 
the ehect of the external classical held on the mentioned physical quantities. The main 
results of the paper are listed in what follows. 

1. Generally, entering the classical held on the interaction together with increasing the 
related parameter may lead to the shift of the maxima amounts of the held entropy 
for two V- and H-type three-level. Increasing the value of 7 leads to the oscillatory 
behaviour of the von Neumann entropy for A-type conhguration. 

2. Maximum values of the DEM between two atoms for H-conhguration are larger than 
those for other two conhgurations. Also, the DEM between two atoms depend on the 
driving classical held. 

3. Since the rank of the reduced density operator of the atomic system containing two 
A-type three-level atoms is three, it is deduced that the PPT is necessary and sufficient 
condition for separability and in this case, the negativity fully captures the entangle¬ 
ment of this system. It turns out that under conditions considered in this work, the 
atom-atom entanglement generated in the A-type conhguration is distillable. 

4. The increase in the external classical held is associated with a slight increment in the 
maximum values of the DEM between diherent subsystems. 
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5. The complete (partial) collapse and revival, as purely quantum mechanical features, 
are observed in Mandel parameter in the absence of the classical field. Also, in the 
presence of classical held the quantized cavity mode is super-Poissonian after certain 
interaction time. 

6 . The numerical results indicate that, no squeezing is seen in y component and squeezing 
occurs in x quadrature in very short time intervals in the beginning of the interaction. 
Also, it is obviously seen that the profundity of squeezing in these regions is decreased 
by the increment of the amplitude of the classical held. 

7. It is illustrated that the amount of considered entanglement criteria as well as quan¬ 
tum statistics and squeezing can be tuned by applying the external classical held 
appropriately via parameter 7 . This, however, is more clear in quantum statistics and 
squeezing as compared with entanglement criteria. 

8 . As a marginal result of the paper, we would like to state that, in the appropriate 
conditions our proposal can be used for the generation of displacement number states. 

At the end of this paper, we mention that this study can be accomplished by considering 
spontaneous emissions and other decoherence processes. We hope to report this work in the 
near future elsewhere. 
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